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Abstract 

We perform a numerical calculation of the total cross section o~(e + e~ — ► 
W+PU - ) as a function of energy, taking into account the finite width of the 
W and the most important radiative corrections. We present these results, 
in tabular form, for several values of M\y. Using these results, we investi- 
gate running strategies for integrated luminosities that might be available at 
LEP 200 and estimate the accuracy to which it will be possible to determine 
the mass of the W by measuring this cross section near threshold. With an 
integrated luminosity of 100 pb" 1 it should be possible to achieve a precision 
of 100 MeV, and with an integrated luminosity of 500 pb -1 , a precision of 
60 MeV. 

PACS numbers: 14.80.Er, 13.38,+c, 13.10.+q 
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I. INTRODUCTION 



The gauge sector of the Standard Model of electroweak physics is described by three 
parameters, which can be chosen to be a, the electromagnetic coupling constant, and Mz 
and M w , the masses of the weak gauge bosons. Of these three, M w is by far the least well 
known. A more practical choice of parameters is a, Mz, and Gf, where Gf = (1.16639 ± 
0.00002) x 10~ 5 GeV" 2 is the Fermi decay constant. The mass of the W, then, can be 
predicted in terms of these three parameters; this prediction depends on some additional 
physics as well, most notably the masses of the top quark and the Higgs boson. Despite this 
dependence, and despite the large uncertainty in the measured value of M w , a comparison 
of Mw to the predicted value is already a stringent test of the Standard Model. 

The W mass is currently measured [|I| to be 80.22±0.26 GeV, where this value is obtained 
by direct reconstruction of W events at pp colliders; it is expected that future experiments 
at the Fermilab Tevatron will be able to reduce this uncertainty to less than 150 MeV. This 
is still, however, a much larger relative uncertainty than in the other fundamental parameters 
of the Standard Model. 

One of the major motivations for upgrading LEP to an energy of ^ ~ 200 GeV is to 
measure Mw more precisely. Several methods have been proposed, including measurement 
of the end point of the leptonic decay spectrum, direct reconstruction of the invariant mass of 
the H^'s decay products H, and measurement of the threshold for W + W~ pair production. 
We will discuss the last of these methods. 

In principle, Mw can be measured simply by determining the end point of the W + W~ 
spectrum; in practice, however, there is no sharp end point, since the W has a finite width, 
with Tw = 2.12 ± 0.11 GeV. Initial-state radiation results in a further smearing of the 
cross section. In the regime \fs ~ 2Mw, which we will call the threshold region, the cross 
section rises steeply as a function of energy, varying from about 2 pb at */s = 160 GeV to 
4 pb at \fs = 162 GeV. A precise measurement of the cross section in this region provides 
essentially the same information as would be obtained from determination of the end point, 
if a sharp end point existed. 

Once the finite W width and initial-state radiation are taken into account, there is 
no simple analytic form for the W + W~ production cross section. Accordingly, we have 
calculated this cross section numerically, and used these results to estimate the precision to 
which it will be possible to measure Mw at LEP 200. 

II. CALCULATION OF THE W PAIR PRODUCTION CROSS SECTION 

At tree level, W + W~ production at an e + e~ collider proceeds by the diagrams in Fig. [I]. 
Evaluation of these diagrams is straightforward; unfortunately, if done naively, it yields an 
unrealistic answer near threshold. The W decays rapidly, and it is incorrect to put it on-shell 
in the final state. 

A more realistic answer Q is given by 
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where do is the cross section for production of two off-shell W's with masses m\ and 777.2, 
and where p is a weight factor, 
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This expression for p differs slightly from that given in Ref. Q], but the difference is numer- 
ically inconsequential. A derivation is given in the Appendix. 

The cross section for pair-production of off-shell W's has been calculated by Muta, 
Najima, and Wakaizumi H, and their result is 
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When both W's are on-shell, i.e., mi = m 2 = Mw, this reduces to the well-known 
form [Efl 
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where (3 is the speed of the W's in the center of mass frame and 
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Numerical evaluation of Eq. [I] is somewhat time-consuming, since it involves two nested 
integrals, both of which are over sharp peaks. Fortunately, it is possible to make a vastly 
simplifying approximation: 
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where the weight factor p is exactly the same as the weight factor p defined in Eq. 
that the W width, IV, is replaced by 2Tw- 

Fig. ^compares the cross section given by Eq. flU) to that given by Eq. (|T0|). The approx- 
imation of Eq. (0), of course, breaks down completely near y/s ^ Mw- This, however, is 
not the domain of interest, and in the region near threshold, y/s ~ 2Mw, the approximation 
is excellent. 

Several other higher-order effects are also significant. The most important are the running 
values of the gauge coupling constants, and initial-state radiation. The first of these can be 
taken into account simply by using the values for the coupling constants renormalized at a 
scale near Mw', we use the values measured at Mz- Inclusion of initial-state radiation 
requires additional work. 

Although initial-state radiation is a purely electromagnetic effect, and is thus suppressed 
by a factor of a, it is nonetheless significant because it is enhanced by a factor of ln(M^/TOg), 
representing the presence of two very different energy scales. Using the formalism of Kuraev 
and Fadin [0, it is possible to sum all orders of initial-state radiation by performing a single 
integral: 
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The second term in the integral represents single-photon hard bremsstrahlung, while the first 
is the result of summing all orders of soft photon emission. This formalism has previously 
been used to include the effects of initial-state radiation to all orders in the calculation 
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of <r(e + e~) — > Z [Q. The difference between this result and the 0{a) calculation can be 
substantial. 

Several groups || have performed full one-loop calculations of a(e + e~ — > W + W~). Our 
calculation is considerably simpler, but contains most of the relevant physics, including 
initial-state radiation beyond leading order in a. 

Fig. [3] shows the results of including initial-state radiation in the calculation of a(e + e~ — > 
Vy -1 "^ - ). Note that initial-state radiation makes a contribution roughly equal in magnitude 
to that of the W's width, and that it has the effect of making the threshold for W pair 
production even less sharp. 

Although there is no longer a sharp threshold, a(e + e~ — > W + W~) still depends strongly 
on the value of the W mass; as shown in Fig. f|, this dependence is strongest for E « Mw- It 
is possible, then, to measure Mw by studying the threshold behavior of a(W + W~). Table | 
presents the same information, for masses within two standard deviations of the current 
central value, in tabular form. 



Note that the behavior of a in the threshold region is mainly due to kinematic effects [10 



so it is proper to use it for a determination of Mw- The behavior of the cross section near 
the peak, at around y/s ~ 220 GeV, is also of interest, but for different reasons. At this 
point, the total cross section depends on delicate gauge cancellations, so it is a sensitive 
probe of the Standard Model gauge structure. 



III. MEASUREMENT OF M w AT LEP 200 



A. Lower bound on the statistical error 



The most favorable situation possible would be if there were no uncertainty in the lumi- 
nosity or the energy of the beam, and if there were no theoretical uncertainties; in this case, 
as seen in Fig. |], measurement of the cross section at even a single point could determine 
Mw, and the error in Mw would be purely statistical. 

The error in the measured value of Mw is given by 
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and, if the error is assumed to be purely statistical, 

5a 



5N 



(14) 



where X = / dC is the integrated luminosity and 5N is the statistical error in the number 
of W + W~ events observed. For the statistical error approaches y/N, and 
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The quantity da/dM can be read off from Table |, but there is a simpler way to obtain 
a rough estimate. Near threshold, as seen in Fig. |], the most important effect of changing 
Mw is simply a shift in the W + W~ spectrum, i.e., 
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where E = y/s/2. Roughly, then, 
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and, measuring the cross section at some particular value of E, the error in Mw is given by 
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Independent of X, then, the optimum energy at which to make this measurement is 
where \da jdE\j ^fa is maximized; this is at E ~ Mw- The precision that may be attained 
is roughly 
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With an integrated luminosity of 1000 pb~ , this would mean 5Mw ~ 30 MeV, which is 
comparable to the precision to which Mz is known. 

As a rough guide to systematic errors, we note that if the measured value of a depends 
multiplicatively on some parameter C, then 
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For example, the measured cross section depends multiplicatively on the detector efficiency, 
and on the luminosity. Certain sources of theoretical error can also, at least approximately, 
be represented this way. Using the values found in Table |, we can rewrite this relation as 
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B. Theoretical error 

There are three potential sources of theoretical error in this measurement: Model depen- 
dence, uncertainty in the input parameters, and incomplete inclusion of radiative corrections. 

In general, a(e + e~ — > lU+iy - ) is model dependent: The calculation we have presented 
assumes the minimal Standard Model, and it will be changed if additional physics beyond 
this framework is included. This effect, however, is most severe far above threshold, where 
delicate gauge cancellations are necessary to preserve unitarity. Near threshold, the behavior 
of the cross section is mainly determined by kinematics, and is less sensitive to the inclusion 
of additional physics. For s/s ~ 160 GeV, model dependence is a negligible source of 
theoretical error. 

Uncertainty in the input parameters is also a negligible source of error. The input param- 
eters used in calculating the W + W~ cross section are the electroweak coupling constants, 
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which, in turn, can be calculated in terms of a, Gf, and Mz- The least well known of these, 
Mz, is still known to better than 0.1%, so this error is also negligible. 

The cross section given in Table [I] is fundamentally a tree-level calculation. We have 
included initial-state radiation, the running of the gauge coupling constants, and the imag- 
inary part of the W's vacuum polarization, but we have neglected several other one-loop 
effects; the dominant such effect is probably the vacuum polarization of the photon and the 
Z. By analogy with physics at the Z peak, we expect that these corrections are less than 
1%. Similarly, we expect that the approximations used in deriving Eq. (|TTD are of this order. 
In the threshold region, the approximation of Eq. ( |TUD is better than 1%, and we neglect it 
compared to other sources of error. 

We will assume, then, that the theoretical error takes the form of an error in the over- 
all normalization, and that it is approximately 1.5%. As discussed in Section [111 A| , this 
corresponds to a systematic error in Mw of approximately 20 MeV. 



C. Systematic error 

One important source of systematic error is the width of the W, IV- In principle, it 
would be possible to determine both Y w and M w simultaneously, in the fit to <j(s). This 
is, however, impractical. As shown in Fig. |5], the shape of the W + W~ cross section is not 
very sensitive to IV, particularly in the region where a(e + e~ —>■ W + W~) is most sensitive 
to Mw- We found, by Monte Carlo simulation, that a simultaneous fit to both Mw and Yw 
resulted in a degraded value of Mw, while failing to provide a more precise value for Yw 
than that which is already known. 

If the Standard Model is assumed to be correct, it is possible to calculate Y w with 
very little theoretical uncertainty. Unfortunately, we cannot use this calculated value in a 
measurement of Mw'- It is, after all, equally possible to calculate Mw in the context of the 
Standard Model, and it would be inconsistent to let Mw vary while using the predicted 
value of Y W - This inconsistency manifests itself even if Y w is computed as a function of 
M W '- Depending on the method of calculation, it scales either as M w or as M w . 

We will, then, simply use the measured value of Y w in our calculation; the uncertainty 
in this value will yield a systematic error in the measurement of Mw- Fortunately, the 
insensitivity of the W + W~ cross section means that this error is small. Performing a nu- 
merical calculation to be described below, and using Y w = 2.12 ± 0.11 GeV, we obtain 
5M W w 20 MeV 

Other sources of systematic error include uncertainty in the luminosity and in the cali- 
bration of the beam energy. From experiences at LEP 100, it is expected that the luminosity 
will be known to better than 1% Jl2] and that the beam energy will be known to within 



20 MeV. Uncertainty in the luminosity affects the overall normalization, while, since we 
are essentially measuring the location of the threshold, an error in the energy calibration 
directly corresponds to an error in Mw 

Finally, performing this measurement requires knowing the efficiency for detecting and 
identifying a W + W~ pair produced nearly at rest. There are two separate issues: The 
probability for misidentifying a W + W~ event as something else, and the probability for 
incorrectly identifying some other type of event as a W + W~ event. Events where at least 
one W decays to a charged lepton and a neutrino are probably distinctive enough that 



7 



misidentification is unlikely, but events where the final state consists of four jets must be 
considered more carefully. 

For our purposes, what is important is not the absolute magnitude of these probabilities, 
but rather the uncertainty with which they are known. This uncertainty, which we denote 
Ae, takes the form of an additional uncertainty in the overall normalization of the measured 
cross section. The parameter Ae can be determined only by detailed detector studies, but 
it is likely that it will be at most a few percent. Depending on its magnitude, this could be 
the dominant source of systematic error. 

We summarize the various sources of systematic and theoretical error in Table [nj. 

D. Realistic estimate of expected precision 

We use a fitting procedure to estimate the statistical and systematic error more precisely: 
We simulate an experimental run by choosing the energies at which measurements will be 
made, and the luminosity that will be devoted to each measurement. Using Poisson statistics, 
we then randomly generate the number of W + W~ events observed at each of these points. 
Finally, we perform a numerical fit of the W + W~ cross section to these randomly generated 
measurements, taking My/ to be a free parameter; an example of such a fit is shown in Fig. 
We can simulate systematic errors by using a different value of Ty^ or of the normalization 
constant when generating the numbers of events than when fitting the cross section. 

The systematic and statistical errors are obtained directly by carrying out this procedure 
many times: The mean of the difference between the fitted value of My/ and the true value 
is the systematic error, while the variance is the statistical error. The different systematic 
errors are added in quadrature. 

We find that, for integrated luminosities between 50 pb~ and 1000 pb~ , and for any 
reasonable assumptions about errors in Ty/ and in the cross section's normalization, the 
optimum strategy is to measure the cross section near threshold, i. e., between 80 and 81 GeV. 
Assuming an integrated luminosity of 100 pb -1 , Fig. [7] shows the statistical error in M w as 
a function of the energy at which the cross section is measured; note that the measurement 
rapidly becomes ineffective as the energy is raised much above threshold. 

The systematic error due to uncertainty in the luminosity and the overall normalization 
of the cross section could be reduced by making additional measurements further above 
threshold, but this systematic error is already less than 50 MeV, which is sufficiently small 
that such an improvement would be negated by the increase in the statistical error that 
would result from lowering the statistics in the region of greatest sensitivity to M w . 

We use the values discussed in Sections |111 B| and |111 (J| to estimate the systematic error. 
Adding the statistical and systematic error in quadrature yields the total error in measuring 
My/, as a function of integrated luminosity; this is plotted in Fig. ||. Until very high lumi- 
nosities are obtained, the measurement is dominated by statistics. Assuming an integrated 
luminosity of 100 pb -1 , we estimate that the total uncertainty will be about 100 MeV. 
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IV. CONCLUSION 



Several methods have been proposed for measuring M w to higher precision than is cur- 
rently available; one of the most promising methods is to measure the threshold dependence 
of the total cross section for W + W~ pair production at an e + e~ collider. 

There is, of course, some tension between this measurement and other physics goals 
of LEP 200: This measurement requires prolonged running in an energy region where the 
W + W~ cross section is only a few picobarns, while many other measurements are best 
made at an energy as close as possible to the peak W + W~ cross section. If it is possible 
to obtain an equally good measurement of M w from reconstruction of hadronic W decays, 
that method may, at least initially, be preferable. 

For these purposes, an elaborate calculation of this cross section is unnecessary: A mod- 
ified tree level calculation includes most of the important physical effects, and is sufficiently 
precise to be compared to the measured cross section for the extraction of M w . We have 
performed such a calculation, and, using these results, we estimate that with an integrated 
luminosity of 100 pb _1 it should be possible to measure M w to a precision of about 100 MeV. 



This work was supported by the Director, Office of Energy Research, Office of High 
Energy and Nuclear Physics, Division of High Energy Physics of the U.S. Department of 
Energy under Contract DE-AC03-76SF00098. 

APPENDIX: PARTICLES OF FINITE WIDTH IN THE FINAL STATE 

If an unstable particle appears in the final state of a process, the naive prescription 
for the cross section is to calculate the cross section for the particle to be produced with 
mass y/s, and then convolve this result with some weighting function p(s) that resembles a 
Breit-Wigner. In this Appendix, we derive a more precise form of this prescription, making 
manifest the approximations that are necessary for the derivation. We will derive a form for 
a final state containing a single W; the generalization to two or more unstable particles is 
trivial. 

The process that actually takes place is the emission of a virtual W, where the final-state 
particles are leptons and hadrons. If box diagrams can be neglected, this process can be 
factorized into the production of a W and its decay, i.e., 



where P M is the part of the amplitude dealing with W production and with any other 
particles in the process, is the W propagator, and J v is the part dealing with W decay. 
Since box diagrams are nonresonant, neglecting them is an excellent approximation. 
Phase space factorizes as well: we can write 
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where qi are the momenta of the final-state particles other than the W, ki are the momenta 
of the W's decay products, p is the momentum of the W, and s = p 2 . The cross section, 
then, may be written 

da = A^W^, (A3) 

where W^ u comes from the W propagator and the Ws decay, and A^ v is everything else. 

By definition, calculating the cross section for production of a W means that we aren't 
interested in the details of what the W decays into, so we sum over all decay channels for the 
W and, for each channel, integrate over the phase space of the decay products. We will deal 
with a single decay channel; the sum over channels presents no additional complications. 
The result, then, is 

= (2tt) 3 J ds A*^A^ j d^ w (h, k n )J* a J p 

= (2tt) 3 / dsA*^A^T a/3 , (A4) 



where 



T 



J d^ w {h, . . . , k n )J^^J v . (A5) 



It is always possible, of course, to write the integrand in Eq. ( |A4] ) as a product, and thus 
to write the cross section in the form 

a = J dsp(s)a (s). (A6) 

This statement, by itself, has no physical content; the real content of Eq. (0) is that it is 
possible to define these functions in such a way that p can be interpreted as a Breit-Wigner 
and <7o can be interpreted as the cross section for production of a W with an unphysical value 
of the mass. It is important to remember, however, that <To is not, in fact, the cross section 
for any physical process, and that its meaning must be defined by explicit construction. 

We will define <To to be the cross section for the production of a W that has an unphysical 
mass but that is still on-shell — that is, that still only has three polarization states. Consider, 
then, the amplitude for the decay of such an on-shell W with momentum p and mass ^/s. 
This W is an ordinary vector particle, so the amplitude is 

M = e;(V~s,p,\)J», (A7) 

where e is the Ws helicity vector, A is the helicity index, and J has the same meaning as 
before. Averaging over initial helicities, 

dT(s) = ^d$ w ■ \-(p"p u - sgnJ\Jv, (A8) 

S O S 

where we have used the fact that 

E^lv^-P- A)e^( v / i,P, A) = -ip^Pu- sg^). (A9) 

A S 
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Performing the phase space integral, 

m = |^ii(pv - «(T)^. (Aio) 

If the masses of all fermions in the final state may be neglected, then the W couples to 
a conserved current, and we may write 

T HV = ^ p v _ S g^)T(s). (All) 

In the case at hand this approximation is permissible: The heaviest fermion that a W can 
decay to is the b quark, and m^/M^ < 0.004. Eq. ( |A10| ) thus simplifies to 

I» = (2vr) 4 v/^. (A12) 

Note, however, that this simplifying approximation is crucial: If the masses of particles in 
the final states cannot be neglected, then no such simple form as Eq. (|2]) exists. 

Similarly, returning to W^", Eq. ( |A4Q simplifies as well. The gauge-dependent piece of 
the W propagator vanishes when contracted with a tensor of the form given in Eq. ([All ), 
so we may write 

A"" = A{s)g^, (A13) 

and, using Eq. ( |A12|) , 

W, u = f ds\A(s)\ 2 ^^Y,^V~s,pA)^,P,\)- (A14) 

The Wfj,u tensor is the only part of the cross section that depends on the details of the W in 
the final state, and this sum over the product of polarization vectors is the form that W^ v 
would have if the final-state particle were a W of mass y/s, so we have established that 

a= I ds\A(s)\ 2 V ^ T ^ a (s). (A15) 

J 7T 

Now consider the form of A(s). At tree level, 

A M = T^T, ( A16 > 

but, including the vacuum polarization of the W it becomes 

— i 

A(S) = s-M^ + 5M^-U ww (sY (A17) 

If M w is taken to be the physical W mass, then Re Hww(Myy) exactly cancels the W mass 
counterterm SM^, so 

A(s) = — ^ , (A18) 

s — M^r — ilm Uww(s) — Re Il ww (s) 
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where II(s) = — n(Mjy). We are only interested in s ~ Myy, and Re fl(s) is a smooth 
function of s that, by construction, vanishes at s = M^- It can thus be absorbed into the 
wave-function renormalization, and we neglect it as a non-leading correction. 
Using the optical theorem, we can show |§ that 

-ImlWOO = VsT(s), (A19) 

where T(s) is the decay rate for an on-shell W of mass ^/s. Finally, then, we can evaluate 
r(s) either by performing the tree- level calculation, or by scaling the measured value. We 
choose the latter method. If the W's decay products are massless, then T(s) scales as y/s, 
so 

r W = v^, (A20) 

where T is the measured width of the physical W. 
Substituting into Eq. (|A15|) , we obtain 



(s - M^f + sr\s) vr 



* = / ds- 2 , — / „ l >- a {s) (A21) 



1 sT 

ds 5 Ti-vois), (A22) 

{s-M^f + s^/M* v kM w 



thus verifying Eq. 
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TABLES 



TABLE I. W + W pair production cross section, including the finite value of and 



the 


effects of initial-state 


radiation, 


for seven 


values of the 


W mass. 


The present value is 


M w 


= 80.22 ± 0.26 GeV. 


The cross 


section is 


a smooth function of Mw, and linear interpo- 


lation should suffice for values of My/ 


in between those presented here. 






y/s (GeV) 






a (pb) 










79.6 GeV 79.8 GeV 


80.0 GeV 


80.2 GeV 


80.4 GeV 


80.6 GeV 


80.8 GeV 


150 


n aai 


fl AO A 


0.402 


0.382 


0.363 


0.346 


0.329 




U.OUo 




0.454 


0.430 


0.408 


0.387 


0.368 


152 


u. ooz 


fl ^A7 
U. 04 I 


0.516 


0.487 


0.461 


0.436 


0.413 


-LOO 


n fi7^ 


fl fi^9 
U.OOZ 


0.593 


0.557 


0.525 


0.495 


0.468 


154 


n 70.^ 
u. < yo 


fl 7*38 
U. 1 OO 


0.689 


0.644 


0.604 


0.567 


0.534 


lOO 


u.yoi 


fl £77 
U.o ( i 


0.812 


0.755 


0.704 


0.657 


0.616 


156 


1 1 fiQ 


1 flfi7 
1 .UO ( 


0.978 


0.901 


0.832 


0.772 


0.719 


1 ^7 

lO 1 


1.4oO 


l.OOO 


1.209 


1.100 


1.007 


0.925 


0.854 


lOO 


i.yoo 


1. ( OO 


1.547 


1.387 


1.252 


1.137 


1.038 


ioy 


z.ooy 


9 

Z.OOO 


2.064 


1.820 


1.615 


1.443 


1.298 


1 fin 

1UU 


O.OOl 


Q 99c; 
0. zZO 


2.831 


2.478 


2.170 


1.908 


1.688 


1 fi1 
1UL 


A 7CIA 
4. I U4 


4. ZOO 


3.812 


3.383 


2.978 


2.609 


2.284 


1 fi9 
1UZ 


^ 79£ 
O. 1 ZO 


^ 90.^ 

o.zyo 


4.853 


4.410 


3.971 


3.540 


3.129 


163 


0.0 ID 


D.zDo 


5.851 


5.427 


4.996 


4.561 


4.128 


1 fiZL 
104 




7 1 f^O 

/ .loy 


6.771 


6.373 


5.968 


5.554 


5.133 


1 fi 1 ^ 
1UO 


O.OZ4 


1 .y 1 u 


7.608 


7.238 


6.860 


6.473 


6.078 


1 fifi 
lOO 


y .UOO 


O. (U/ 


8.369 


8.024 


7.672 


7.311 


6.943 


1 fi7 


y .ooy 


Q ^7£ 

y.o/o 


9.061 


8.744 


8.409 


8.073 


7.730 


168 


1U.Z04 


Q QQ9 

y .yyz 


9.694 


9.391 


9.082 


8.767 


8.446 


1 fiQ 

iuy 


lU.OOU 


1U.OOZ 


10.272 


9.987 


9.696 


9.400 


9.099 


170 


1 l.OOU 


1 1 flfift 
1 l.UDo 


10.803 


10.533 


10.257 


9.979 


9.696 


1 71 

-L I 1 


1 1 78Q 

ii. i oy 


1 1 . 041 


11.289 


11.034 


10.774 


10.510 


10.241 


172 


12.212 


11.976 


11.736 


11.494 


11.247 


10.997 


10.744 


173 


12.601 


12.376 


12.148 


11.917 


11.683 


11.445 


11.204 


174 


12.961 


12.746 


12.528 


12.307 


12.084 


11.858 


11.629 


175 


13.294 


13.088 


1 9 878 
1Z.O 1 


1 9 fifiS 
IZ.OOo 


19 /IK/1 
IZ.^404 


12.238 


12.020 


176 


13.601 


13.403 


13.203 


13.001 


12.797 


12.590 


12.381 


177 


13.885 


13.695 


13.503 


13.309 


13.113 


12.915 


12.715 


178 


14.148 


13.965 


13.781 


13.595 


13.406 


13.216 


13.024 


179 


14.391 


14.215 


14.037 


13.858 


13.678 


13.495 


13.310 


180 


14.616 


14.446 


14.275 


14.102 


13.928 


13.752 


13.575 


181 


14.824 


14.660 


14.495 


14.328 


14.160 


13.991 


13.820 


182 


15.016 


14.857 


14.698 


14.537 


14.375 


14.212 


14.047 


183 


15.193 


15.040 


14.886 


14.730 


14.574 


14.416 


14.257 


184 


15.357 


15.209 


15.060 


14.909 


14.758 


14.605 


14.451 


185 


15.508 


15.365 


15.220 


15.075 


14.928 


14.780 


14.631 


186 


15.648 


15.508 


15.368 


15.227 


15.085 


14.942 


14.798 


187 


15.776 


15.641 


15.505 


15.368 


15.230 


15.091 


14.952 


188 


15.894 


15.767 


15.630 


15.497 


15.364 


15.229 


15.093 


189 


16.002 


15.875 


15.746 


15.621 


15.487 


15.356 


15.225 


190 


16.102 


15.977 


15.852 


15.727 


15.600 


15.477 


15.345 
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TABLE II. Sources of systematic error in the measurement of Mw by study of the threshold, 



as discussed in Sections IIIC and HID; errors are added in quadrature. Many of these errors take 
the form of uncertainty in the overall normalization, which allows us to use Eq. ©. We are unable 
to estimate the uncertainty in detection efficiency, so we parameterize it as Ae. In our plots, we 
assume Ae = 2%. All of these errors are sufficiently small that the measurement will be dominated 
by statistical error. 



Effect 


Effect on normalization 


Error in Mw 


Coupling constants 


negligible 


negligible 


Model dependence 


negligible 


negligible 


One-loop corrections 


1% 


17 MeV 


Initial-state radiation 


1% 


17 MeV 


Smearing approximation 


negligible 


negligible 


Total theoretical error 


1.5% 


20 MeV 


T w 




20 MeV 


Luminosity 


1% 


17 MeV 


Beam energy 




20 MeV 


Detection efficiency 


Ae 


(17 MeV) x Ae 


Total systematic error 




45 MeV 
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FIGURES 



FIG. 1. Tree-level Feynman diagrams for W + W pair production at an e + e collider. Ex- 
change of a virtual Higgs boson also contributes, but is negligible. 

FIG. 2. Comparison of an exact tree-level calculation of a(e + e~ — > W + W~), in which both 
W's in the final state are integrated over Breit-Wigners, to an approximate calculation, in which 
only one W is integrated over a Breit-Wigner, but using a width twice the physical value. Note 
that this is an excellent approximation except in a domain so far below threshold that the cross 
section is unobservably small. The peak at 45 GeV is due to the decay of a Z to two virtual W's. 

FIG. 3. Calculation of a(e + e~ — > W + W~). The solid curve is the tree-level cross section 
where both final-state W's are on-shell, the dashed curve is the cross section where the W's are 
allowed to be off-shell, and the dotted curve includes initial-state radiation as well as the finite W 
width. Note that the finite width of the W and the initial-state radiation are effects of comparable 
magnitude. 

FIG. 4. Calculation of <r(e + e _ — > W + W~), including the finite width of the W and the 
effect of initial-state radiation. The solid curve is for M\y = 80.0 GeV, the dashed curve is for 
Mw = 80.2 GeV, and the dotted curve is for M\y = 80.4 GeV. The dependence on M\y is strongest 
for E rs M w . 

FIG. 5. Calculation of a(e + e~ — > W + W~), including the finite width of the W and the 
effect of initial-state radiation. The three curves represent different values for T\y: the solid curve 
is with F\y = 1.9 GeV, the dashed curve is with T\v = 2.1 GeV, and the dotted curve is with 
T\v = 2.3 GeV. The measured value is = 2.12 ± 0.11 GeV. Note that for values of the width 
close to the measured value, the form of the cross section is not very sensitive to T^- 

FIG. 6. Simulation of an experimental determination of M\y, using an integrated luminosity 
of 100pb -1 . Measurements are made at 80.1 GeV, 80.5 GeV, and 80.9 GeV. At each point, the 
number of observed events is randomly chosen, using a Poisson distribution whose mean is the 
cross section times the integrated luminosity. The points, with 1-cr error bars, represent these 
simulated measurements of a at those three energies, and the solid curve, which corresponds to 
M w = 80.31 GeV, is the best fit. 

FIG. 7. Statistical error in the measurement of Myy, in MeV, as a function of the energy at 
which a(e + e~ — > W + W~) is measured. Note that the statistical error is quite high unless the 
measurement is made at the W + W~ threshold. This plot is generated assuming an integrated 
luminosity of 100 pb -1 , but this qualitative feature is true independent of luminosity. The total 
cross section at this energy is roughly 0.3 pb, compared to a peak value of more than 16 pb. 
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FIG. 8. Estimated error in the measurement of Myy, in MeV, as a function of integrated 
luminosity. The solid curve is the total error, and the dotted curve is the statistical error. Note 
that until very high luminosities are obtained, the measurement is dominated by statistical error. 
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Figure 1 
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